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O ! Abstract 

, The major obstacle to a supersymmetric theory on the lattice is the failure of the 

Leibniz rule. We analyze this issue by using the Wess-Zumino model and a general 
Ginsparg- Wilson operator, which is local and free of species doublers. We point out 
that the Leibniz rule could be maintained on the lattice if the generic momentum 
k^ carried by any field variable satisfies \ak^\ < 5 in the limit a — > for arbitrarily 
small but finite 5. This condition is expected to be satisfied generally if the theory is 
finite perturbatively, provided that discretization does not induce further symmetry 
\ breaking. We thus first render the continuum Wess-Zumino model finite by applying 

' the higher derivative regularization which preserves supersymmetry. We then put 

this theory on the lattice, which preserves supersymmetry except for a breaking in 
interaction terms by the failure of the Leibniz rule. By this way, we define a lattice 
{Sj ' Wess-Zumino model which maintains the basic properties such as U(l) x U(1)r 

symmetry and holomorphicity. We show that this model reproduces continuum 
theory in the limit a — > up to any finite order in perturbation theory; in this sense 
all the supersymmetry breaking terms induced by the failure of the Leibniz rule are 
irrelevant. We then suggest that this discretization may work to define a low energy 



Oh: 



. effective theory in a non-perturbative way. 

> : 

^ 1 Introduction 

There are basically two different motivations for defining a field theory on the lattice. 
The first is to regularize a divergent theory and simultaneously define the theory in a 
non-perturbative sense. The second is to define a discretized version of a theory, which 
is finite in continuum perturbation theory, so that one can apply the numerical and 
other techniques in a non-perturbative way. Though this second aspect is not commonly 
discussed in the context of lattice theory, we want to show that this second aspect may 
be essential in putting supersymmetric theories on the lattice 1 . 

There are several difficulties to define supersymmetry on the lattice. The most notable 
and difficult issue is the failure of the Leibniz rule[f2j. To be explicit, we have on the lattice 

-{f{x + a)g(x + a) - f(x)g(x)) 
a 

1 If one applies a discretization to superstring theory,for example, it also corresponds to a discretization 
of a perturbatively finite theory. 
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= -(/0 + a) - f(x))g(x) + f{x)-{g(x + a) - g(x)) 
a a 

+ J-(/(x + a) - f(x))^(g(x + a) - g(x)) (1.1) 

namely the "lattice version of the Leibniz rule" is given byf\ 

{V{fg)){x) = (Vf){x)g(x) + f{x)(Vg){x) + a(V f){x){Vg)(x). (1.2) 



This shows that the breaking of supersymmetry by the lattice artifact is formally of 
order 0(a), but actually the breaking is of order 0(1) if the momentum carried by the 
field variables is of order 0(1 /a). To recover the conventional Leibniz rule, a necessary 
condition for the momentum variable is 



\ak^\ < 5 (1.3) 

for a — > with arbitrarily small but finite S. Here is a generic momentum carried by 
any field variable in the Feynman diagrams so that the last term in the lattice Leibniz 
rule (1.2) is neglected to give 

(V(fg))(x) = (Vf)(x)g(x) + f(x)(SJ g)(x). (1.4) 

This requirement is expected to be satisfied if the theory in continuum is finite in a 
perturbative sense so that all the momentum variables in Feynman diagrams are finite 
and thus infinitesimally small measured by the lattice unit 1/a in the limit a — > 0, provided 
that the lattice discretization does not introduce further symmetry breaking terms. 

If the above condition is satisfied, all the supersymmetry breaking terms for finite 
lattice spacing, which are induced by the failure of the Leibniz rule, are expected to 
be irrelevant in the sense that those supersymmetry breaking terms vanish in the limit 
a — ► 0. In the context of supersymmetric Yang-Mills theory on the lattice, an argument 
to the effect that all the supersymmetry breaking terms are irrelevant was given in the 
past|| , though the N = 1 supersymmetric Yang-Mills theory is not finite and thus the 
basic reasoning is completely different. In the context of the Wess-Zumino model^j, we 
would like to show that a sensible lattice discretization, which is based on the presently 
available technique, may be to first render the continuum Wess-Zumino model finite by 
applying the higher derivative regularization. This higher derivative regularization is 
known to preserve supersymmetry in a perturbative sense§. We then apply the lattice 
discretization to this regularized continuum theory. 

In this paper, we present a detailed analysis of the above procedure for the Wess- 
Zumino model in the framework of perturbation theory, and then suggest that this scheme 
may work in a non-perturbatice sense also. We utilize the Ginsparg- Wilson fermion 

2 If one uses the symmetric difference, 

(V(fg))(x) = (f(x + a)g(x + a) - f(x - a)g(x - o))/(2o) 

= (Vf)(x)g(x) + f(x)(Vg)(x) + a(Vf)(x)[(g(x + a) - g(x))/a] + a[(f(x) - f(x - a))/a](Vg)(x) 
and still the limit a —> is not smoothly defined in general. 
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operators which are local and free of species doublers flj - []I0| . In the course of this analysis, 
we clarify some of the subtleties appearing in the Ginsparg- Wilson operators when utilized 
in the present context. 



2 Wess-Zumino model on the lattice 

The Wess-Zumino model is the simplest supersymmetric model in 4-dimensional space 
time, and the non-renormalization theorem was first discovered in this model: If renor- 
malized at vanishing momenta, all the potential terms including mass terms do not receive 
any (even finite) renormalization, except for a uniform wave function renormalization, up 



to all orders in perturbation theory [I I 



As for the previous studies of the Wess-Zumino model on the lattice, see, for example, 

I3-IH. 



2.1 Lattice Lagrangian 

We define the Wess-Zumino model on the lattice byp~5[ 



£ = l -X T C^ 1 ,Dx+ l -mx T Cx + 9X T C(P + <PP + + P-(t> ] P-)x 

l 1 5 L 



-$D*D<j> + Fi-jF + m[F<p + (F0) f ] + g[F<p 2 + (F0 2 ) f ] 
\x T C^r- X + \rn X T C X + g X T C(P + <pP + + P-^P-)x 

-(f>^(f> + F^F + m[F<p + (F0)t] + g{F<f) 2 + {F^}. (2.1) 
a 1 5 



Here the hermitian operator 

H = a^D = H ] (2.2) 
satisfies the general Ginsparg- Wilson relation 

75 # + #75 = 2H 2k+2 (2.3) 

with a non-negative integer k, which implies 

l5 H 2 = ( l5 H + Hj 5 )H - H( lb H + H l5 ) + H 2 l5 = # 2 75 , (2-4) 

and 

T 5 = 75 - H 2k+1 (2.5) 

which satisfies HT 5 + T 5 H = 0. An explicit form of H is given in Appendix. When we 
have H 2 and T 2 in the bosonic terms, we adopt the convention to discard the unit Dirac 
matrix in H (see Appendix). Note that Y\ = I — H 4k+2 , and H 2 = 1 just on top of the 
would-be species doublers. 
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The projection operators are defined by [10 

P± = ~(l±7 5 ), P±= l -{1±%) (2.6) 

with 75 = 75 — 2H 2k+1 which satisfies t| = 1. 

A salient feature of our Lagrangian is that it is invariant under the continuum chiral 
transformation, except for the mass term, if one performs simultaneously a suitable phase 
rotation of the fields <fi and <jy . We here note the relation which follows from the defining 
relation of H (see also Appendix) 

T 5 H/a = -( l5 H-H^ k ) = ±[ l5 ,H}<xY S -^. (2-7) 
a la a 

We then have 

{ 7 5,r 5 tf} = i{75,[75,tf]} = (2.8) 
which suggests that the fermion kinetic operator satisfies 

{75, ^H} = {75, ^T 5 H} = (2.9) 
1 5 1 5 

by using [75, T§] = 0. The factor T 5 in the fermion kinetic termP] (1/T 5 )H vanishes at 
the momentum corresponding to the would-be species doublers [[L5[] []17|] , but this non- 
locality is compensated for in the present supersymmetric theory by the corresponding 
singularity in the term F^F. Note that variables {x, 4>, F} are treated as components of 
a single superfield in supersymmetric theory. In fact one can confirm that the partition 
function of the free part of the Lagrangian gives unity and thus the factor is cancelled 
among the component fields: 

J VxV4>V^VFVF^exp{ J [l x T C±-^ X + \m X T C X 
_0t^0 + F \1_ F + m[F + (F0)t]]} 

d 1 c; 
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v x V(t>v^VF'v{F')UM f i-x T c- — x + -™x T c x 

J Z 1 5 d Z 

-0t^ _ t (mr5) 2 + (F / )t J_^ ]} 

a 1 5 



det [j^f + m\ Jdet [±] det [f + mT 



det[f + (mr 5 ) 2 ] det[^] det[f + K 5 ) 2 ] det^ 
{det[^]det[(f) 2 + (mr 5 ) 2 ]} 



1/4 



det[5 + (mr 5 ) 2 ] det[^j 



(2.10) 



3 This structure of the fermion kinetic term is required to define the Euclidean Majorana fermion|T( 
in a consistent manner in the presence of chiral symmetric Yukawa couplings [la]. 
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if one recalls that both of r§ and (— ) 2 are proportional to a 4 x 4 unit matrix and that 
this unit matrix is neglected in the bosonic sector appearing in the denominator. 

We shall also confirm that perturbation theory is well defined without any singularity, 
though the fermion propagator vanishes at the momenta corresponding to the would-be 
species doublers. In the non-perturbative formulation, the factor I/T5 may be compen- 
sated for by rescaling the variables F and F\ as will be shown later. 

Our convention of the charge conjugation matrix is 

C^C- 1 = -(7 M ) T , (2.11) 

C-feCT 1 = 7?, (2.12) 
tfC =1, C T = -C (2.13) 

and the Ginsparg- Wilson operator satisfies 

c^r.c- 1 = ( 75 r 5 ) T , 

CDC- 1 = D T . (2.14) 



2.2 Supersymmetry 

If one defines the real components by 



±={A + iB), F^-^(F-zG) (2.15) 



the above Lagrangian is written as 



C = X T C--H X -— 2 [AH 2 A + BH 2 B] 
I 1 5 a 2a z 



+ l [F fl F + G fl G] 
1 



+ ^mx T Cx + m[FA + GB] 

-^9X T C(A + z l5 B) X + -^9[F(A 2 - B 2 ) + 2G(AB)}. (2.16) 

The free part of the action formed from this Lagrangian is confirmed to be invariant under 
the "lattice supersymmetry" transformation^ 

5 X = -T 5 -H(A - t l5 B)e — (F — z 75 G)e, 
a 

5 A = e T C X = X T Ce, 

SB = -ie T C l5X = ~iX T C l5 e 

5F = e T CT 5 -Hx ~ x T CT 5 -He 
a a 

5G = te T CT 5 -H l5X ~ -*X T CT 5 -# 75 e (2.17) 
a a 



4 We assume that our kinetic operator T^H satisfies the proper charge conjugation symmetry 
CT^HC^ 1 = (T^H) 7 " , which includes an operation corresponding to partial integration. 
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with a constant Majorana-type Grassmann parameter e. Note that the order of the 
operators is important in these expressions. The second expressions of 5F and SG need 
to be treated carefully when these variations are multiplied with other field variables 5 . If 
one recalls the correspondence to continuum theory in the naive limit a — > 

-H <-> 75 V> 
a 

r 5 ^ 75 , (2.18) 

the above transformation defines a lattice generalization of the continuum supersymmetry 
transformation || || . 

For example, the variation of the kinetic terms under the above transformation is given 

by 



5 / C 



kin 



= I { X T C^-H[-T 5 -H(A - t l5 B)t — (F — i 7s G)e] 

J 1 5 CL d 

-- 2 AH\ T C X ] - - 2 BH^e T C l5X ] 
cr a 1 

+F^[e T CT 5 -H X ] + G^[ie T CT 5 -H^ X ]} 

1 g CL 1 g CL 

= [ { X T C- 2 H\A - i l5 B)e - X T C^-H{F - i l5 G)e 
J a z 1 5 a 

-- X T CR\A - t l5 B)e + /C<-^r 5 -tf (F - i^G)t} = (2.19) 
cr 1 5 a 

which is consistent with (2.10). Here we used (2.8). 

The lattice supersymmetry variation of the interaction terms is given by 



^ J C-int 

I { 72 9XTC{A + i -y* B )[- r *l H ( A ~ l ^B)t ~{F- i l5 G)e] 
+l=g{FA[e T C X } - FB[-ze T C l5X } + 2GA[-ze T C l5X } + 2GB[e T C X ]} 

+^=g{(A 2 - B 2 )[e T CT 5 -H x ] + 2(AB)[ie T CT 5 -H l5 x)}} 
V2 a a 

J {l= gx T C(A + i 7s B)[-r s ±H{A - i l5 B)e] 

+^=g{(A 2 - B 2 )[e T CT 5 -H x ] + 2(AB)[ie T CT 5 -H %X }}} 

J {l=g X T C(A + tl5 B)[-T 5 ^H(A - z l5 B)e] 



'For example, in the presence of a scalar field A(x), we have J Ae 7 CT$-Hx = J X T CF5^H(Ae) 
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+^=g{[ X T CT 5 -H(A 2 - B 2 )e] - 2[i X T CT 5 -H l5 (AB)e]}} 

-^=g X r C(AT 5 -HA - BT 5 -HB)e 

+ ^-g X T C[AY,-HB + BT 5 -HA] l5 e 
V2 a a 



/< 



+^=g{[ X T CT 5 -H(A 2 - B 2 )e] - 2i[ X T CT 5 -H(AB)j 5 e}}}. (2.20) 

Here we used the relation 7 5 r 5 ii = — T 5 Hj 5 (2.8). If the operator T 5 H/a (2.7) satisfies 
the Leibniz rule, the above variation of the interaction terms vanishes. We thus encounter 
the notorious issue related to the Leibniz rule, which is basically the lattice artifact. 
The propagators for perturbative calculations are given by 



a 2 



H 2 + (amT 



5 J 



\2 



rr2p2 

(FFt) = (-) 



if 2 + (amT 



5 2 



H 2 + (amf 



5y 



\2 



<*«x T W = H^^rr. = ( - )75 r 5¥T ^- 75 (2.2i) 

and other propagators vanish. When we have if 2 and r 2 in the bosonic propagators, we 
adopt the convention to discard the unit Dirac matrix in H 2 . Note that T\ = 1 — ff 4fc+2 . 
If one uses the identities 

p + <P(x)P + = p + 0(x)p +75 r 5 , 

P-<f>\x)P- = P^(x)P^ l5 T 5 , (2.22) 

the Feynman rules in the present scheme are essentially identical to those in the previous 
calculation The one-loop level non-renormalization theorem when remormalized at 
vanishing momenta is thus satisfied, though the kinetic terms receive non-uniform finite 
renormalization in addtion to uniform logarithmic renormalization[|T^]. 

The holomorphic properties in a naive sense are preserved in our Lagrangian[]. As 
for the U(l) x Ur(1) charges, where Ur(1) stands for the R-symmetry, we first write the 
potential part of the Lagrangian as 

Cpot = l -m{P +X ) T CP + x + mF<P + x T CP + g<pP + x J rgF4> 2 (2.23) 

+ l -m\P_ X ) T CP-X + (rnFtp) 1 + X T CP-(g^P- X + Q?F0 2 ) f 

6 The decisive factor to ensure the non-renormalization theorem is supersymmetry, while U(l) x Z7r(1) 
symmetry and holomorphicity provide additional constraints. 
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and we assign [19 



0= (1,1), 
^ = (1,-1), 
P + X = i= (] 
m = (-2,0), 

9 = (-3,-1) 



1,0), 



(2.24) 



by regarding m and g as complex parameters. Here £ is a two component spinor in the 
representation where 75 is diagonal. For m = g = 0, our Lagrangian preserves these 
charges if one recalls 



and P-X = m the representation where 75 is diagonal. 

3 Higher derivative regularization on the lattice 

In the perturbative treatment of the above lattice Lagrangian (2.1), the higher order 
diagrams generally break supersymmetry because all the momentum regions contribute 
to the loop diagrams; it is not easy to preserve supersymmetry in this strict sense ( i.e., 
for all the momentum regions ) on the lattice because of the failure of the Leibniz rule. 

A way to resolve this difficulty may be to apply a higher derivative regularization on 
the lattice. By this way, one can transfer all the divergences to the infrared divergences 
measured by the lattice unit 1/a. In the infrared region, the momenta in loop diagrams are 
constrained to the momentum region in continuum theory, for which the lattice artifact 
such as the failure of the Leibniz rule could be negligible in the limit a — > 0. 

The higher derivative regularization in the present lattice Lagrangian is implemented 




HP 



(2.25) 



by 



C = \x T C^D 

l 1 5 



(D^D + M 2 ) 
M 2 



X + C 



{D^D + M 2 ) 
M 2 



X 




X 



H 2 (H 2 + (aM) 2 ) , 1 (H 2 + (aM) 2 ) 

a 2 (aM) 2 r§ (aM) 2 



(3.1) 




+g X T C(P + <pP + + P-^P-) X + g[F<P 2 + (F0 2 )t]. 
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The U(l) x Ur(1) symmetry and holomorphicity are preserved in this regularization. 
The propagators for perturbative calculations are given by 



> 

(FF+) 



t a 2 (aMf 



H 2 + (amT 5 ) 2 H 2 + (aM) 2 

H 2 T 2 (aM) 2 
~~ C > H 2 + (amT 5 ) 2 H 2 + (aM) 2 

a 2 mT 2 (aM)' 



<F0> = <Ftyt) = (_) 



H 2 + (amT 5 ) 2 H 2 + (aM) 



2 



T a (aM) 2 

( X (y) X (x)C) = H75r 5 g + Qmr5g2 + (aM)2 75 

aT 5 H - a 2 m(r 5 ) 2 (aM) 2 
" H 2 + (amT 5 ) 2 H 2 + (aM) 2 



(3.2) 



and other propagators vanish. When we have H 2 and T 2 in the bosonic sector, we adopt 
the convention to discard the unit Dirac matrix in H 2 . Note that T5H + HY§ = and 
T 2 = 1 — H 4k+2 . Here M is a new mass scale which may be chosen to be 

1 > (aM) 2 > (am) 2 . (3.3) 

One can confirm that the free part of this Lagrangian with higher derivative regular- 
ization (3.1) is still invariant under the lattice supersymmetry transformation (2.17) by 
noting [rs, H 2 ] = 0, while the interaction terms are not modified by the higher derivative 
regularization. 

3.1 One- loop tadpole and self-energy corrections 

It has been shown previously that the superpotential is not renormalized in the one-loop 



level even for a finite a when renormalized at vanishing moment a [|18|j. This conclusion still 
holds in the present model with higher derivative regularization. One can also confirm 
that the cancellation of tadpole diagrams is still maintained even for a finite a in the 
one-loop level in the present modelQ. 

It is instructive to analyze the tadpole diagrams in some detail. The scalar tadpole 
contribution to the fields <fi and is given by 

2s[(F0)0+(Ft0t)0t] 

a 2 mr 2 (aM) 2 f a 2 mr 2 (aM) 2 



-29 Ji 

-2mga 2 (<j) + ^) J 



H 2 + (amT 5 ) 2 H 2 + (aM) 2 H 2 + (amT 5 ) 2 5 H 2 + (aM) 2i 
H 2 + (amT 5 ) 2 H 2 + (aM) 2 



d A k T 2 (k) M 2 

-2mg($ + ) / ^ ( 2vr )4 #2 (jt ) + ( am r 5 (k)) 2 H 2 (k) + (aM) 2 ( } 



intuitively, this one-loop cancellation arises from the fact that the interaction terms, when one of 
4> or (p^ is set to a constant, are reduced to the effective mass terms which are invariant under lattice 
supersymmetry. 



9 



where we chose the basic Brillouin zone at 



— < K < - (3-5) 
a a 

and rescaled the integration variable as ak^ — > k^. If one considers the limit a — > in 
the above integral, one obtains a logarithmic divergence from the region k^ ~ 0. But it 
is important to recognize that the entire region of the basic Brillouin zone gives a finite 
contribution in the above integral even at the limit a — > 0. This is the peculiar feature of 
the one-loop tadpole diagrams in the present minimal higher derivative regular iz at ion|, 
and all other diagrams receive non-vanishing contributions only from the infrared region 
k^ ~ in the limit a — > 0. In any case, it is confirmed that the above scalar tadpole 
contribution is precisely cancelled by a fermion tadpole contribution even for a finite a. 

Since the renormalization of kinetic terms was not uniform without the higher deriva- 
tive regular ization|TS[, we here analyze the kinetic terms and associated quadratic and 
logarithmic divergences in the one-loop level in more detail. The one-loop correction to 
the "kinetic term" FF^ is given by 

||[F(0) 2 + (F(0) 2 )t] 2 



^[4F<#t)^t)jrt] 



2g 2 a*F[ J 



(aMf 



x 



H 2 + (amT 5 ) 2 H 2 + (aM) 2 
1 (aM) 2 



H 2 + (amT 5 ) 2 H 2 + (aM)' 



]F f (3.6) 



which is logarithmically divergent. 

The integral for FF* is written in more detail as 



2g 2 



d 4 k 1 (aM) 2 

(2vr) 4 H 2 (k + ap) + (amT 5 ) 2 (k + ap) H 2 (k + ap) + (aM) 2 



1 (aM) 2 
X H 2 (k) + (amT 5 ) 2 (k) H 2 (k) + (aM) 2 ^ ^ > 

where we rescaled the integration variable as ak^ k^ by choosing the basic Brillouin 
zone as in (3.5). In this integral, if one chooses the integration domain outside the infrared 
region 

5 > k^ > — 5 for all ji (3.8) 

for arbitrarily small but finite S, the integral vanishes for a — > since we have no infrared 
divergences^. In this analysis, the absence of species doubling in the Ginsparg- Wilson 

8 If one considers the higher derivative regularization with the factor ((H 2 + (aM) 2 ) / (aM) 2 ) 2 instead 
of (H 2 + (aM) 2 ) I (aM) 2 in (3.1), even the one- loop tadpole diagrams receive non- vanishing contributions 
only from the infrared region in the limit a — » 0. This stronger regularization may be necessary for a 
numerical simulation. 

9 This and following analyses are extended to 8 = (aM) e with sufficiently small positive e. 
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operator is essential: Namely, H 2 ~ 1 for the momentum domain of the would-be species 
doublers. 

The non-vanishing contribution to the above integral in the limit a — > is thus given 

by 



* d 4 k 1 (aM) 2 

s (2tt) 4 H 2 (k + ap) + (amT 5 ) 2 (k + ap) H 2 (k + ap) + (aM)' 



1 ( flM ) 2 1 (39 ) 

H 2 (k) + (amY b ) 2 (k) H 2 (k) + (aMf 1 1 ' ; 

for arbitrarily small but finite 5. We can thus write this integral as 

r* d 4 k 1 (aM) 2 



8 (2tt) 4 (k + ap) 2 + (am) 2 (k + ap) 2 + (aM) 2 
1 (aM) 2 . 



F + (am) 2 A; 2 + (aM) 2J 



s/a (2tt) 4 (jfe + p) 2 + m 2 (k + p) 2 + M 2 
1 M 2 , 



X 



/c 2 + m 2 k 2 + M 2i 

2r /■<» d A k 1 M 2 



oo (2vr) 4 (k + p) 2 + m 2 {k + p) 2 + M 2 
1 M 2 

X k 2 + m 2 k 2 + M 2 ^ (3 ' 10) 

in the limit a — > 0. Here we used H 2 (k) ~ A; 2 for |fc| < 5 by recalling the rescaling 
a/c M — >• fc^. See Appendix. This last expression is identical to the continuum result in the 
higher derivative regularization. 

The crucial aspect of this analysis is that the momentum variables are constrained 
to the infrared region in the limit a — > 0. Namely, the typical momentum variable is 
constrained to be 

\ak»\<5 (3.11) 

for arbitrarily small but finite 5, which is a necessary condition for the validity of the 
Leibniz rule on the lattice. 

The correction to the kinetic term of the scalar particle (f> by the fermion loop diagram 
is given by 

2 2 r ,f (aM) 2 T 5 H (aM) 2 

9 a rq> H 2 + (mar 5 ) 2 H 2 + (aM) 2 H 2 + (mar 5 ) 2 H 2 + (aM) 2 1 } 

where the symbol Tr includes the integral over the loop momentum as well as the trace 
over Dirac matrices. The quadratic divergence and the logarithmic divergence associated 
with the mass term in this expression, when evaluated at vanishing external momentum, 
are cancelled by the scalar loop diagram (3.14) given below even for a finite a. 
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By analyzing the infrared structure in the limit a — > 0, one can confirm that the 
expression (3.12) is reduced to the continuum expression in the higher derivative regular- 
ization 

H f d ' k ) ^+ 1 M2 ,v ^ I M " m<n 

9 J (27r) 4 (k + p) 2 + m 2 (k + p) 2 + M 2 k 2 + m 2 k 2 + M 2 ' 1 ' 

The one-loop self-energy of the scalar particle produced by the scalar particle loop 
is given by 

||[F0 2 + (F0 2 )t][F0 2 + (F0 2 )t] 

-» |[8#t(^) W t)] 

-H 2 Y 2 (aM) 2 a 2 (aM) 2 



9 f -\ 

= 21 y [800 ^ 2 + ( am r 5 )2 #2 + ( qM )2 #2 + ( am r 5 )2 #2 + ( flM ) 2 J 

= -4g a 00ty r 5H2 + {aM)2H2 + {amT5)2H2 + {aM)2 (3.14) 



2! J H 2 + (amT 5 ) 2 H 2 + (aM) 2 /J 2 + (amT 5 ) 2 H 2 + (aM) 2 ' 

The first term is quadratically divergent, and the remaining term is logarithmically di- 
vergent. These terms, when evaluated at vanishing external momentum, precisely cancel 
the corresponding fermion contributions (3.12) even for a finite a. 
The fermion self-energy correction is given by 

gay +[J + ^ am p 5 ^ 2 H 2 + ( a M) 2 H 2 + (amr 5 ) 2 H 2 + (aM) 2 * 

2 37 P f r 5# {aM) 2 1 (aM) 2 

+ 9 a V - [ J H 2 + (amT 5 ) 2 H 2 + (aM) 2 H 2 + (amr 5 ) 2 H 2 + (flM) 2 ' +t 

(3.15) 

where we used the relation 

P ± T 5 H = T 5 HP T . (3.16) 
By analyzing the infrared structure in the limit a — > 0, we again have 

27 P r y ^ ^ i m 2 

^ +[ i (27r) 4 (A; + p) 2 + m 2 (A; + p) 2 + M 2 A; 2 + m 2 fc 2 + M 2j 



, , p , d 4 A; ^ M 2 1 M 2 

+ 5 ^ (2tt) 4 (fc + p) 2 + m 2 (k + p) 2 + M 2 k 2 + m 2 k 2 + M sJ 



(3.17) 



We have to examine if a universal wave function renormalization is sufficient to remove 
the divergences from these 3 contributions (3.10), (3.13) and (3.17) in the limit a — > 0. 
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One can in fact show that a uniform subtraction of logarithmic infinity ( for M — > large) 
renders all these expressions finite when renormalized at vanishing momentum. For the 
fermion contribution to the scalar kinetic term (3.13), we first rewrite it as 



Ag 2 J — - 1 </. t (p)0(p)3! j dad/3d-fd55(l - a- /3- 7 -5)x 

-p 2 {l -a-f3)M A 



(3.18) 



[p 2 (a + (5){l-a- (5) + a{k 2 + m?) + (5{k 2 + M 2 ) + 1 (k 2 + m 2 ) + 5(k 2 + M 2 )} 4 



and similarly for the fermion self-energy correction (3.17). We then renormalize all the 
kinetic terms at p = by using the relation 



In the one-loop level, we can thus maintain supersymmetry including renormlization fac- 
tors 10 in the limit a — > 0. In other words, all the supersymmetry breaking terms for finite 
lattice spacing should vanish for a — > 0. 

3.2 Two and higher-loop diagrams 
3.2.1 Tadpole diagrams 

We start with the analysis of tadpole diagrams. One can confirm that the tadpole di- 
agrams for the auxiliary field F in the two-loop level precisely cancel even for a finite 
lattice spacing a, and similarly for the field F\ 

For the scalar field 0, we have 4 tadpole diagrams in the two-loop level. These diagrams 
do not quite cancel for a finite a due to the failure of the Leibniz rule for general momenta. 
But for non- vanishing contributions in the limit a — >• 0, one can reduce the Feynman 
amplitudes to those of the continuum theory with higher derivative regularization, which 
are then shown to cancel precisely. We can thus maintain the vanishing tadpole diagrams 
in the two-loop level. 

For example, we have a two-loop tadpole contribution arising from a fermion loop 
diagram with a scalar exchange correction, which contains a quadratic divergence in a 
naive sense, 




(3.19) 



[a(k 2 + m 2 ) + [3(k 2 + M 2 ) + 7 (P + m 2 ) + 5(k 2 + M 2 )] 4 ' 




10 



Here we are repeating the known analysis in continuum 
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a 2 m(r 5 (p)) 2 M 2 

X H 2 (p) + (amT 5 (p)) 2 H 2 (p) + (aM) 2 

aV 5 H(k+p) M 2 
X H 2 (k + p) + (amT 5 (k + p)) 2 H 2 (k + p) + (aM) 2 
a 2 M 2 
# 2 (fc) + (amr 5 (fc)) 2 # 2 (fc) + (aM) 21 l ' ; 

where we used the rescaled variables ap^ — > p M and afc^ — > fc^. By analyzing the limit 
a — > 0, one can confirm that only the infrared regions 

Ippl < 8, \kfil < S for all \i (3-21) 

with arbitrarily small but finite 8 give a non- vanishing contribution 

W l & I 5 —tr{ aV " H{T>) ^ 



8 (2tt) 4 Js (2tt) 4 l H 2 (p) + (amT 5 (p)) 2 H 2 (p) + (aM) 2 
a 2 m(T b (p)) 2 M 2 



x 



H 2 (p) + (amT 5 (p)) 2 H 2 (p) + (aM) 2 

aT 5 H(k + p) M 2 

H 2 (k + p) + (amT 5 (k + p)) 2 H 2 (k +p) + (aM) 



2 



a 2 M 2 

X H 2 (k) + (amT 5 (k)) 2 H 2 (k) + (aM) 2 >' (3 ' 22) 

By rescaling the momentum variables to original ones and considering the limit a — » 0, 
one recovers the continuum result in the higher derivative regularization. Other tadpole 
amplitudes in the two-loop level are analyzed similarly. 



3.2.2 General diagrams in two and higher-loop level 

Similarly, one can confirm that the non-vanishing parts of all the Feynman amplitudes 
in the two-loop level in the limit a — > are reduced to those of the continuum theory 
with higher derivative regularization. Namely, all the loop momenta in the Feynman 
amplitudes are constrained to be in the infrared region, where the Leibniz rule is satisfied. 

To the extent that the non-renormalization and other good properties are maintained 
in the continuum theory with higher derivative regularization ]|], our lattice regularization 
thus reproduces all the good properties of the supersymmetric Wess-Zumino model in the 
limit a — ► 0. 

One can extend this analysis up to any finite order in perturbation theory, since the 
power counting in this regularized theory is effectively super convergent; the higher order 
diagrams are thus more ultra-violet convergent and thus less sensitive to the lattice cut-off 
for a — > 0. All the supersymmetry breaking terms induced by the failure of the Leibniz 
rule are thus shown to vanish in the limit a — > 0. 
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3.3 Non-perturbative treatment 

As for the non-perturbative treatment of our model, one may first perform the path 
integral over the Majorana fermion which produces the Pfaffian 



r 1 H(H 2 + (aM) 2 ) (H 2 + (aM) 2 ) . , . . 

2/IJ2 („A/r\2\ 1 /rr2 i /„ » / f\2\ 



r r t li H 2 (H 2 + (aM) 2 ) , , 1 



(# 2 + (aM) 2 



M) 2 r T 2 (aM) 2 
m ( J ff 2 + (aMf), /n (ff 2 + (aM) 2 ) A+1 . ~ ,„, 2 , +111 



where in the second expression we rescaled the field variables as 

F -> TF, F f -> F f r (3.24) 



with T = \JT 2 = \/Y — H Ak+2 . In the last expression in (3.23) we have no singularity 
associated with l/r 5 . This path integral (or after performing the path integral over F 
and Ft ) may be evaluated non-perturbatively. Since the Wess-Zumino model is not 
asymptotically free, the non-perturbative result in the continuum limit a —>■ may be 
defined with a finite M, which provides a finite mass scale to specify the renormalized 
parameters; one can thus prevent the coupling constant from increasing indefinitely. In 
this sense our possible non-perturbative formulation, which is inferred from perturbative 
considerations, is consistent. 

If the above path integral (3.23), when evaluated non-perturbatively, gives a well- 
defined result in the limit a — > 0, it is expected that the path integral defines a theory 
which incorporates the quantum effects up to the energy scale M. The path integral may 
then be effective in defining a non-perturbative low energy effective theory^ for 

|?V|~m«M (3.25) 

where p M is the typical external momentum carried by field variables. 

n If one should be able to evaluate the continuum theory with higher derivative regularization in a 
non-perturbative way, one would obtain the same result. 
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4 Discussion 



We have discussed a way to ensure the Leibniz rule for the supersymmetric Wess-Zumino 
model on the lattice. The basic observation is that the lattice Leibniz rule is reduced 
to that of continuum theory if the generic momentum carried by any field variable 
is constrained in the infrared region \ak^\ < 5 for arbitrarily small but finite S in the 
limit a — > 0. A way to ensure this momentum condition is to apply the higher derivative 
regularization to the Wess-Zumino model so that the theory becomes finite up to any 
finite order in perturbation theory. On the basis of the analysis of Feynman amplitudes, 
we have shown that this is in fact realized in our lattice formulation which incorporates a 
lattice version of higher derivative regularization. All the supersymmetry breaking terms 
induced by the failure of the Leibniz rule thus become irrelevant in the sense that they 
all vanish in the limit a —>■ 0. We suggested that this mechanism may work even in a 
non-perturbative sense. 

In this analysis, it is crucial that the Ginsparg- Wilson operator is free of species 
doublers so that the lattice operator H is of order 0(1) in the momentum regions of the 
would-be species doublers 12 . Besides, the Ginsparg- Wilson operator maintains some of the 
basic symmetries such as chiral symmetry, U(l) x Ur(1) symmetry and holomorphicity. 

In conclusion, we have presented a possible way to maintain the Leibniz rule on the 
lattice. It is yet to be seen if this analysis is extended to supersymmetric Yang-Mills 
theories on the lattice p0|-||21||, though we naively expect that supersymmetric theories 
which are perturbatively finite ( such as N = 4 theory ) may be put on the lattice 
consistently. One may also want to find a more drastic way to overcome the difficulty 
associated with the Leibniz rule, which might lead to a completely new understanding of 
lattice regularization. 



A Ginsparg- Wilson operators 

An explicit form of the general Ginsparg- Wilson operator H 0, which satisfies the algebra 
(2.3), is given in momentum space by 

| fc + l I fc+1 I fc + 1 I k Sf 

H(a Pll ) = lb (-)^(^=)—i{^H^ + M k )—^-^Hl l -M k )—H^} 
1 ^ H w a 



7s(^) ^ (-7=) + M k) »+i - (^F\k] - M k ) mh-i i $ 



(A.l) 

where k is a non-negative integer and 

F W = (s 2 ) 2k+1 + Ml 

12 It is possible to implement the present mechanism for the model based on the Wilson fermion|l2|], if 
one defines the higher derivative regulator suitably by including the Wilson term. The non-perturbative 
analysis would, however, become more involved since the basic symmetries such as chiral symmetry and 
holomorphicity are spoiled. 
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with 



M k = E(i 

A 1 



(A.2) 



s M = sin ap M 
c M = cos ap^ 

^= 7 M sin ap M . (A. 3) 



This operator is known to be local and free of species doublers|22|], and this operator for 
k = is reduced to Neuberger's overlap operator f8j. Our Euclidean Dirac matrices are 
hermitian, (7**)* = 7 M , and the inner product is defined to be s 2 > 0. Note that H 2 ( 
and consequently Tf = l — H 4k+2 ) is independent of Dirac matrices. The parameter mo 
is constrained by < mo < 2 to avoid species doublers, and 2mQ fc+1 = 1 gives a proper 
normalization of H, namely, for an infinitesimal p M , i.e., for \ap^\ <C 1, 

H ~ - T5 ai i>{\ + O(ap) 2 ) + 75 ( 75 az ^) 2fe+2 (A.4) 

to be consistent with H = j^aD; the last term in the right-hand side is the leading term 
of chiral symmetry breaking terms. 
We thus have 

H 2 = 1 (A.5) 
just on top of the would-be species doublers, for example, ap^ = (n, 0, 0, 0), and 

H 2 ~ (aprf (A.6) 

for \apfj,\ 1 independently of the parameter k. 
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